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o\ : 

We reexamine in detail the problem of fermions interacting with a non-Abelian random vector 
, potential. Without resorting to the replica or supersymmetry approaches, we show that in the 

hmit of infinite disorder strength the theory possesses an exact solution which takes the form of a 
logarithmic conformal field theory. We show that the proper treatment of the locality conditions in 
the SU(2) theory leads to the termination of the multifractal spectrum, or in other words to the 
\^ ' termination of the infinite hierarchies of negative-dimensional operators that were thought to occur. 

Based on arguments of logarithmic degeneracies, we conjecture that such a termination mechanism 
should be present for general SU{N). Moreover, our results lead to the conclusion that the previous 
replica solution of this problem yields incorrect results. 
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There now exists a great deal of evidence that the wavefunctions in disordered systems near a locahzation- 
delocalization transition exhibit multifractal bchaviom'. In other words, the moments of the density of states scale 
not according to one given fractal dimension like in a simple fractal, but rather according to an infinite set of scaling 
exponents. Since the most interesting features of multifractality are beyond the reach of perturbation theory, one 
''O • needs to employ non-perturbative approaches. In this context the exactly solvable problem of Dirac fermions in a 
Ph random vector potential in two dimensions has proved to be a very instructive training ground. This model was first 
^ , considered in [Q in connection with the integer quantum Hall transition. Recently exact results for the full multifractal 
spectrum for Dirac fermions interacting with a random magnetic field (i.e. an Abelian vector potential) have been 
obtained from a Gaussian field theory in an ultrametric space |^], and reproduced by mapping to a random energy 
^ ; model |]. 

ly-j , The non-Abelian random vector potential case, originally introduced in [Q as an effective model for low-energy 
' excitations around the nodes of a d-wave superconductor, has very recently appeared in physical problems of a 
very different nature, namely in non-Hermitian quantum mechanics There, it appears as the effective theory 
of a system driven by a strong imaginary vector potential in the limit of weak disorder. Non-Hermitian quantum 
mechanics describes such problems as anomalous diffusion in a random media |^ and statistical mechanics of flux 
lines in superconductors ||7[]. The multifractality is observed in impurity-averaged correlators of local moments of 
^ wavefunctions, for which expressions for the non-Abelian part need to be computed. Although these correlators were 

investigated in [^,^, we wish to report some new and surprising aspects of the problem that were overlooked. An 
, earlier presentation of some initial results can be found in ||l0[| . 
I • The non-Abelian Random Vector Potential (RVP) has already been treated by both standard approaches for 
^ , disordered systems, namely the replica and SUSY ones. The SUSY treatment can be found in ||], while the replica 
Q ' treatment was put forward in and carried through in It was first recognized in this publication that the 
^ , associated conformal theory should be of a special class, namely a logarithmic one. However, both of these lines of 
ILj^ ■ investigation suffered from various setbacks. First of all, the SUSY approach seemingly generated an infinite series 
of operators with negative conformal dimensions for SU{N) randomness. On the other hand, the replica approach 
gave a different physical content, coming from the fact that the conformal blocks of the theory were slightly different. 
There was consequently an urgent need to provide a better comprehension of this system by performing some more 
careful and extensive investigations. 

As was pointed out in pd| ], the disorder averaging in the non-Abelian RVP problem can be performed without 
using the SUSY or replica approaches in the limit of vanishing frequency and infinite disorder strength. This limit 
corresponds to the conformal limit of the theory, which was studied in ||] . The RVP model thus is a critical disordered 
system, which makes it theoretically very interesting in its own right. In addition, the logarithmic nature of the CFT 
involved calls for better understanding. 

Usually, CFTs are synonymous with power-law dependence of physical correlators. Much work has been done 
about the unitary minimal models, for which there exists a unitary finite-dimensional representation of the Virasoro 
algebra. However, recently it has become clear that logarithmic dependence can appear for models outside of this 



1 



class. The first instance in which logarithms were introduced in conformal correlators is to be found in jlj]. It was 
first recognized in JT^ that this type of logarithmic behaviour was associated to non-diagonalizability of the Virasoro 
operator. Subsequent applications of logarithmic operators were found in critical disordered models Cp^i and 

Cp^q models [|l5|-p^, critical polymers and percolation |16-lJ], gravitationally dressed CFTs (l9j, two-dimensional 



magnetohydrodynamic turbulence pOlM , D-brane recoil in string theory |22-G4] , and were recently shown to corre- 



spond to a novel bulk excitation in the Quantum Hall state 1 25 . Further new developments about logarithmic CFTs 
can be found for example in [2^-|3l|] . 

This paper aims to resolve some important remaining issues related to fermions interacting with a non-Abelian 
random vector potential in the conformal limit. Our main result is to show that the correct treatment of the full 
conformal field theory that is obtained from the SUSY or direct treatment (which coincide) leads to the termination 
of the parabolic multifractal spectrum before the scaling exponents r* (g) reach their maximal value. They were 
previously thought to follow this parabolic law to arbitrarily negative values. 

The paper is organized as follows. In section II we introduce the model and present the conformal limit. In section 
III we manipulate the theory at the critical point and recall how to average over disorder without SUSY or replicas. In 
section IV we explicitly construct the free-field formulation of the theory for = 2, with which multipoint correlators 
are to be calculated explicitly in section V. Section VI discusses the fusion rules and the operator algebra of the 
logarithmic conformal field theory, which are shown to have great impact on the scaling behaviour to be expected. 
The treatment for > 2 is presented in section VII together with a discussion about the extension of the N — 2 
results to this domain. In section VIII we discuss the consequences of our results on the multifractal spectrum. Section 
IX discusses the consequences on the replica solution, and we finish with conclusions in section X. 

II. THE MODEL 

We consider N species of Dirac fermions living in a 2+1-dimensional space and interacting through a disordered 
vector potential belonging to an su{N) algebra A, to which they are coupled minimally. The disorder allows for 
hopping between the different species. Since the vector potential is time- independent, different Matsubara frequencies 
do not couple, and can be treated independently by a Euclidean two-dimensional theory with explicit frequency 
dependence. In fact, for a given realization of the disorder, the partition function is given by the fermionic path 
integral 

Z[A^] = yp[*, ^] e-5[*.".'4^1 (1) 

with the Dirac action 

S[^,Lu,A^]= [d^x^{x)[I(» ^-iuj + i4]^{x) (2) 



(since we are in a two-dimensional Euclidean space, we take the Pauli matrices as Dirac 7 matrices, i.e. jjL 
Disorder-dependent single-particle Green's functions can then be represented by the path integral 



G(x, ^; A,) = ^ jv[^, ^ vI/(a;)*(y)e-^[*'--^-l (3) 

Physical quantities are obtained by performing the disorder averaging procedure on products of Green's functions. 
We use the distribution functional 

\nP[A^]^~^ j d?xTvA^{x)A^{x) (4) 

representing the usual ^-correlated Gaussian white noise for the random vector potential. 
The impurity-averaged Green's function then reads 

G{x,y,u:)^G(x,y-Lo-A^)^ [ VA,Gix,y;cj; A^)P[A,] ^ -i ( VA,^^ ^{x)^{y)e-^^'''^^^-^ (5) 



'Z[A,,j 

In the limits of infinite disorder strength g ^ 00 and of vanishing frequency a; — > 0, the theory becomes conformally 
invariant. We will formulate and solve for correlation functions at this conformal point, and later use renormalization 
group arguments to infer the scaling behaviour of physical quantities like the density of states and its local moments 
away from criticality. We will use the notation G{xi, ...,x„) to denote the impurity-averaged n-point correlator at 
zero frequency. 
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III. THE THEORY AT THE CRITICAL POINT 



In the following sections, we will be concerned with the step-by-step formulation of the theory at the conformally 
invariant critical point, i.e., for infinite disorder strength and vanishing frequency. The principles of this approach 
have been formulated by Bernard and Mudry et. al. ^ . In this work we show that the path integral is perfectly 
well-defined and negative-dimensional operators do not appear in the physical operator OPEs. 

We start by separating the fermionic action into chiral parts: 

S[^,A^] = S+[^,A+] + S-[^,A^] 

S±[^,A±] = J d^x^i{x)[d±+iA±ix)]^±{x) (6) 

where we have used the holomorphic and antiholomorphic derivatives and fields {2d = 9_, 29 = 9+) 

d± ^ di ± 182 

A±=Ai± 1A2 (7) 

where now A± g su'~'{N) — sl{N,C), the complex extension A'^ of A. This translates, at the level of the path 
integration over the vector fields, into transforming the double integral of Ai over su{N) into a single integral over 
su''^{N). The non-compact nature of the SU'~^{N) group manifold, in contrast to SU{N), will be of crucial importance 
later on. 

We parametrize the vector fields by fields g± belonging to the complex extension G"^ of the group G — SU{N): 

A±{x) ^ id±g±{x)gZ\x) (8) 

The reality condition A^_^_{x) — A^{x) translates into g\_{x) = gZ^{x). From now on, we will use the notation 
9+{x) = g{x). 

Let us study, for a little while, the problem at fixed disorder. We can observe that the transformations 

*±(x) ^ g±{x)^'±{x), ^iix) ^ ^'J{x)g^\x) (9) 

completely decouple the fermions from the random vector potential. It maps the minimally coupled Dirac action for 
^ fermions into the free Dirac action for ^E"' fermions, i.e. 

S[^,A^,] = j d?x[^'\d+-^\ + ■^'^_d-'^'_] (10) 

generating the correlators 



The Jacobian for this transformation. 



{^'_{-Z,)^'l{-Z2)) ^ (11) 

27rzi2 



Det[^ 

has the very important property of being proportional to the partition function at fixed disorder (up to the anomalous 
determinant Det[^, which is irrelevant for computation of correlation functions but contributes to the total central 
charge), thus cancelling it when computing the correlations for fixed disorder, i.e., 

G{xu ...) - jvA.^J V[^, i']vl/...e-^[*^^l = / / ^^I*'. S^'-e-^^*'! (13) 

This removes the need to invoke either the replica or supersymmetry methods to perform explicitly the disorder 
averaging. This also implies that we are not dealing with a free energy functional, and that disconnected diagrams 
do appear in a perturbation expansion. 
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The original disorder averaging procedure involved two su{N) algebra integrations (for the two components of the 
vector field). We transform these two integrations into an integration over a "gauge" sector and a "coset" sector, 
by writing A± — ^(f'^^^e) _|_ ^(coset) rpj^^ change of variables (^) transforms these two integrals in a single group 
integration over the complex extension SU'-^{N), which we can again separate into "gauge" and "coset" integrations 
by writing g± = uH±,u £ SU{N),H± € SU^{N)/SU{N) (the reality condition states then that H^^ = ij|). This 
procedure induces a non trivial Jacobian in the path integral: 

VAf, = VA1VA2 = X)A(sa«ge)p^(coset) ^ jyj^{gauge)jyjj^2NW[H^H] ^--^^^ 

where VH is the integration measure over the coset SU'~'' (N) / SU (N) , and N corresponds to the dual Coxeter number 
of SU{N). W[h] is the WZNW functional [dB = S) for the field h = W H on level k = -2N (note that the 

topological term then has positive sign): 

W[h] = ^ / cPxTiid^hd^h'^] + fd'x e,,,xTr[d^hh-^d,hh-^dxhh~'] (15) 

87r Jg IZTT Jg 

Thus, the partition function for N species of Dirac fermions interacting with an su{N) random vector potential has 
completely factorized in the product of two independent sectors: a free fermion (disorder independent) part, and an 
SU'~^{N) part in which all disorder dependence is concentrated. 



IV. COSET FACTORIZATION AND FREE FIELD FORMULATION OF THE WZNW MODEL 

The important point to notice from the previous formulas is that, since the WZNW model depends only on the 
combination H which is invariant under left-multiplication by w € SU{N), i.e. H uH, it is advantageous 
to proceed to the factorization of SU'^'{N) into the subgroup SU{N) and the coset SU'^' {N)/ SU{N). This can be 
explicitly performed with the help of the Iwasawa decomposition of SU'-^{N) |^^. The Haar measure on SU^{N) 
then factorizes into the product of the Haar measures on SU'^ (N) / SU{N) and SU{N). 

For simplicity, we specialize now to = 2. Our derivation follows essentially the approach to WZNW models as 
free field theories to be found in |3^. For g e SU^{2), the Iwasawa decomposition takes the form 











[ 









(16) 

with u £ SU{2) and G R, G C. The integral measure for this parametrization is then given by (see derivation in 
the Appendix) 

dg = du e'^'^d(f>dfj,^dfi- (17) 

du being the Haar measure on SU{2). 

The field appearing in the Jacobian (p^), being SU{2) independent, can be expressed in terms of the fields ((/), fi-) 

as 

Using this parametrization in the Jacobian (|lj) gives 

kW[h] = ^ / d^x[id^(t)f + e'^'f'd-^i+d+iJ,-] (19) 
in J 

as the action for the SU'~^ [2) / SU [2) fields (0, /i+,/i_). For negative fc, this action is positive definite. Once again: 
even though the Jacobian (|l^) induces a WZNW model with a negative level, which is not well-defined at the level 
of the path integral for a group manifold with a positive metric, the fact that the WZNW field h belongs to the 
coset SU'~^ (2) / SU {2) which has a negative-definite metric makes the coset action a positive-definite functional. This 
is reminiscent of the coset constructions on non-compact manifolds to be found in p5[ | , where the coset formed of 
the non-compact group and its maximally compact subgroup forms a well-defined unitary conformal theory when one 
considers a negative level for the underlying Kac-Moody algebra. 
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The WZNW model has the property of being invariant with respect to Kac-Moody chiral current algebras. The 
Kac-Moody currents, given by (for negative level) 

J{z) = ~kdhh-\ J{z) = -kh-^dh (20) 

have, in terms of the new fields, the representation 

J{z) = J+(T- + J-CJ+ + Ha^ 

(21) 

J- = -ke^'''dfi+ 

H = -k[d(t)- ^l-e^'l'^^x+\ 

(22) 

H ^H^ J± = J=F^ (23) 

In the path integration over (0, /i^, /i-), the presence of the factor e^*^ in the measure as well as the presence of the 
combination e^'^dji^ in the expression for the current J{z) make the change of variables 

0J^{z) = -ke^'l'^^l+ (24) 

very convenient, since it brings us to a path integration over free fields. 

This and the properly regularized anomalous determinant associated to this transformation |32] 



Det[e2^9] ^ {Bei[e-'''^de'*d]Y'^ (25) 



which shifts the action for (p by 



1 

47r 



d^x[2{d^,(l)Y + TZcj)] (26) 

171 J 

modify the coset action to 

S'[0,cj_,^_] = -5- / (fx[{-k-2)[dAf+n<l)] + ^ I d^xu^d+n^ (27) 
An J 2tt J 

i.e. the anomaly shifts k ^ k + 2 = —q^ (recall that in the given case k = —4, = 2) in front of <j) and introduces the 
Riemann curvature TZ of the manifold in the action for (jj, which will change its conformal properties. By rescaling 
(f) — > -;^4>' , the action and Operator Product Expansions for our free fields read 

oTT J q Ztt J 

(28) 

(/)'(zi)0'(z2) = - In2;i2, 0'(zi)(^'(z2) = - lnzi2, ^_(zi)w_(z2) = — (29) 

2^12 

in which we have performed the separation of 0' in terms of holomorphic and anti-holomorphic parts (j)'[z,z) = 
(f>{z) + 4>{z), valid at the level of correlation functions. 

The proper expression for the currents in terms of our free fields is then given by (^3|), but only after performing 
the same transformations as above on (f>, i.e. the shift kcj) {k + 2)(j) and the rescaling (p ^ '^2q^' ' 

J-{z)=u;^{z) 

J+{z) = V2q^i-{z)d<j>'{z) ~ kdfi^{z) - filizP^iz) 

H{z) = ^qdp'{z) - ^l-{z)Lo^{z) (30) 

In all of these, the usual normal-ordering procedure is implied. The currents have the correct OPEs for generators in 
the Cartan-Weyl basis of SU{2)k, i.e. 
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i7(zi)J±(z2) = — J±(Z2) + ... 
2:12 

J+(Z1)J-(Z2) - + —H{Z2) + ... 

2^12 Z12 

Hiz,)H{z2) = ^ + ... (31) 

Thus, what we have done is the foUowing: by using the Iwasawa decomposition (Jl^), we have exphcitly shown that 
the level k = —2N = —4 cosct SU'^ {2}/ SU{2) possesses a Wakimoto free field representation with Cartan-Weyl 
generators carrying an affine SU (2) current algebra with level analytically continued to negative value. 
The stress-energy tensor 

nz) = ^[lij+j- + j~,n+H^] (32) 
2 



is, in terms of our free fields. 



T{z) = -l-d(l)'d(l)' - ^d^(f>' -uj^dfi^ (33) 
2 V2(7 



The central charge is readily checked to be c = 3k/{k + 2) as for the normal case of positive k (in our case, c = 6). 
The 0' part of T{z) is that of a Dotsenko-Fateev free field with imaginary background charge ap ~ */(2'z). Vertex 
operators Va{z) = e°"^ '■^^ then have holomorphic conformal weight = —(a/2) (a + \/2/q). 

Primary operators Vp can be found by requiring the correct OPEs with T{z) and the generators of the Kac-Moody 
algebra (|3^). A highest- weight finite-dimensional representation is available, composed of the 2j + 1 operators Vjm{z) 
given by 

Vjm{z) - e4^^'(^)[-/i-(z)]^+'" (34) 
for which the OPEs with the currents have the correct SU (2) form 

Tfi 

H{zi)Vj„,{z2) - —Vj„,{z2) (35) 

2^12 

These primary operators all possess identical conformal weights 

a,V2 jjj + l) jjj + 1) 
h,„, - - 2 (a + — ) = ^ = -jTfY (36) 

(note that /i_ is a weight zero operator, as can be seen from its OPE with T{z), and consequently does not contribute 
to the weight of Vjm). Weights are thus negative in any spin-j representation for our value k — —4. 

To calculate correlation functions, we follow the general procedure outlined in This starts with the introduction 
of dual operators Vjm , generalizing the usual dual vertex operators of the Coulomb-gas formalism. The highest-weight 
operator Vj-j is given by 

VS-_,(z)=e-f(^+^>'(z)..^+^^-(.) (37) 

Acting on this with J+ produces the other states in the representation, but there is a slight setback: they do not 
have such a simple form as (37). To obtain a correlator involving another dual state than (|3^), the trick is to use the 
conformal Ward identities and derive it from another correlator involving (|3^). For future reference, we here give the 
two dual operators in the fundamental spin- 1/2 representation: 



vi(z)^yi/2._i/2 = er*^'^V(z) 



V+{z) = Vi/2.i/2^ '^V2qd(t>\z)ujl{z)- ijL_{z)ujt{z)-l2duo-{z)ujl.{z) e v^^"^ (38) 

The general prescription for the calculation of correlators then requires that we replace one of the fields by its 
dual, and use screening charges if necessary to satisfy the neutrality conditions, which are now modified to 
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= --^=-3, N^_-N^_^3 (39) 

with a understood as the numerical coefHcient of the vertex operator, used as in ev^?'^ , and where iV(^_ and 
count these respective operators in the correlator. Note that this last condition is equivalent to the requirement that 
the correlator be a singlet state, e.g. that it contains an equal number of + and — indices in the spin-1/2 case. 

The last ingredient we will need is the screening operator. It can be determined in two ways: first, it must commute 
with all current generators up to a total differential; but also, it can be found from the change of variables ([24|). 
Inverting this, we can see that over every noncontractable cycle on the Riemann surface, we must have 

dC^_(C)e~^'^'('^) = -k(fdn+ = (40) 



Enforcing this condition can be done by putting in the path integral the (infinite) product, over all closed cycles, of 
5 functions 

Y[S{<j^ uj-e-'f^') ^ J d\exp{iX ^ iv.e-'r't'') (41) 

Expanding in powers of A then produces the screening charge insertions, which we will denote by the operator Q{z). 

We have explicitly constructed the holomorphic operators of the theory. An identical procedure can be done for 
the antiholomorphic part, using instead of ( p^ ) the conjugate transformation 

uj+{z) = -ke^'I'd^i- (42) 

This way, primary highest-weight states of the antiholomorphic currents can be constructed, as well as their duals, 
and a conjugate screening charge can be found. We will denote all of these by V , etc. 

The matrix field h{z, z) can then be expressed, in correlation functions, in terms of the tensor product of left and 
right primary spin-1/2 fields as 

h{z -z) - ( ^-^'^^--^'^ (43) 



V+{z)V-{z) V+{z)V+{z) 
The two-point functions of primary fields 

(T4(^i)^t(^2))=4' (44) 

give for the coset two-point correlator 

{haibdzi:Zl)h^^\^{z2,Z2)) = (Jaias'^biba 1^12 1^^^ (45) 

We have thus seen how coset correlations can be simply obtained by calculating the relevant conformal blocks of 
our primary fields. This will be done for the four-point function in the following section. 

V. MULTI POINT CORRELATION FUNCTIONS 

Let us start by moving back to our original physical problem, and see what we are now capable of doing. By 
considering the limiting case g —^ oo, we have decomposed our original fermions into free fermions ip' together with 
SU'~^ (N) / SU{N) coset operators h. Thus, by considering correlators of the local operators 

0-\z,z) = *t ^vl,_^^ KIK,'^-, (46) 



we explicitly sum over the SU{N) indices, leaving only the coset operators once the decoupling transformation (|^) 
has been done. The SU (N) path integration then simply factorizes out of the effective generating functional. 
The two-point function 
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{0{l)0-\2)) ^ ^(Tr[Ml)/.-i(2)]) ^ (47) 

shows that the operator O has conformal weights (1/8, 1/8) (i.e. a spinless operator of dimension 1/4). This reproduces 
the previously known resuh 

Let us now turn to the four-point function. By performing the fermionic contractions, we get 

(0(1)0-1(2)0(3)0-1(4)) ^ ^^(Tr[/^(l)/i-l(2)]Tr[/^(3)/l-l(4)]) + 

+ ^- {Tr[h{\)h-\2)h{i)h-\A)]) +(2^4) (48) 

By using (^), we can see that all correlators appearing in (|8|) can be expressed in terms of the three vertex correlators 

Ci = {V+{z^)VAz2)V+{z^)V^{z^)) 
C2 = {V+{zi)V+{z2)VAz:i)V^{zi)) 
C3 = {V-{zi)V+{z2)V+{z3)V-{z4)) 

These can be calculated using the prescription of the last section: we replace the operator in position 4 by its dual 
(note that, in the definitions above, we have performed an overall spin flip, which does not change the correlators but 
allows to use the more simple expression V-), and insert the necessary screening charges (here, only one). We use the 
projection invariance to perform the analytic mapping (zi, Z2, 23, Z4) (0, z, 1, 00), z = f^^f^, for which we get 

a(zi...Z4) = [zi3Z24]^/*C,(z), aiz)^ lim zU^a{0,z,l,Zoo) (49) 

Performing the integrals (with Ci{z) — [z{l — z)]'^/^Ci{z)) yields 

Cr(z)=^^(3/2,3/2;2;z) 

(7i^(z)=F(3/2,3/2;2;l-z) 

C2"(z)=^^(l/2,3/2;l;z) 

C2'(^) = ^^^(l/2,3/2;2;l-z) 
C^{z)^^Fil/2,S/2;2;z) 

6-3^(2) =F(l/2, 3/2; l;l-z) (50) 

in which the two independent solutions for each correlator are labeled by the integration contours a : (0, z), 5 : (z, 1) 
used in the screening. The b contour solutions have logarithmic behaviour near z — 0, whereas the a ones have 
logarithmic behaviour near z = 1. We are thus in the presence of a logarithmic conformal field theory, whose first 
examples appeared in | |l^ , [T^ . 

The "full" correlator, in which the z dependence has been included, has to be built by combining the conformal 
blocks is such a way as to satisfy the crossing symmetry constraints, as well as the requirement of single- valuedness. 
Crossing symmetry is easily satisfied in view of the extremely simple monodromy properties of the solutions ([s^). 
The constraint of single valuedness (for example, around z = 0) forces the cross-multiplication of a and b solutions 
respectively in the z and z sectors, without which terms of the form Inzlnz would appear, which obviously are to be 
excluded in a single- valued correlator. 

In fact, by explicitly multiplying out ( p3[ ) and using the above solutions, we find, using the standard identities 

^^(3/2, 3/2; 2; z) + ^^^(1/2, 3/2; 2; z) = ^(1 - z)F(3/2, 5/2; 2; z) 
i^(3/2, 3/2; 2; z) - F(l/2, 3/2; 1; z) = ^zF(3/2, 5/2; 3; z) 

(51) 

that the four-point correlator can be written as a function of the conformal blocks 
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F^{z) = 2V2{1 - z)F{3/2, 5/2; 2; z) 
F^{z) = (1 - z)F{3/2, 5/2; 3; 1 - z) 
F^{z) = zi^(3/2,5/2;3;z) 

/'l(z) = 2^/2zF(3/2,5/2;2;l-z) (52) 

These can be verified to be the conformal blocks that one would obtain by solving straightforwardly the SU{2)^4 
Knizhnik-Zamolodchikov equations. We have obtained them using a somewhat different route, through the coset 
parametrization. Thus, even though naively the SU{2)^4 WZNW model seems ill-defined at the level of the path 
integral (since the kinetic term, for a negative level, has the wrong sign !), the fact that we are in reality working 
on the coset SU^{2)/SU{2) which has been shown above to carry the same current algebra makes the use of the 
KZ equations valid. Even though we have proved this in detail only for N = 2, we can expect that the same will 
be true for iV > 2 as well. One very important remark to make at this point is that these conformal blocks do 
not reproduce the conformal blocks obtained within the replica approach |^ . Although the operator dimensions are 
correctly reproduced, the correlators and OPEs are incorrect. We will return to this point later on. 
The full coset correlator 

Ha^,,,b^{zi,Z2,Z3,Z4) = \zi3Z24:\^^^Ha^,,,b^{z,z) 

= (/laibi(2;i,Zi)/l-^^^(z2,Z2)/la3b3(2:3,Z3)/lfr/a4(^4,^4)) (53) 

can now be built as mentioned above, by solving the monodromy problem, and insuring single-valuedness in the 
complex plane. The correlator can be projected as usual onto the singlets 

^1 — ^aia2^a3a4 ^2 ^aia4^a2a3 

Hai...b4{z,z) ^ ^ lJjHij{z,z) (54) 

with, by single-valuedness, 

i/y(z,z) = a\z{l ~ z)\^/^[F^^iz)F^{z) + (a ^ b)] (55) 

where a is some constant chosen to satisfy consistency with the normalization of the two-point function. 
The full four-point correlator for the local operators (^), with free fermion contribution, finally reads 



(0(l)...0-i(4)). 





3/2 1 


2122^342142^23 


213224 





+ - - } Hai,,,b^(Z,Z) (56) 



212234 214Z23 J I 212234 214223 



where P and P are projectors onto the singlets, i.e. P{Ii)P{Ij)Ha^,,,bi — Hij- 

For further discussion we shall need the four-point correlator of the physical field M (related to the local density 
of states), which is composed of the sum of the two operators (|4£ 



M{z,z) =Tv^^ = 0{z,z) + 0-^{z,z) (57) 

The correlator we are interested in is made up of the different permutations of (^) that appear when we make use of 
(^). In fact, it reads 

{M{l)M{2)M{'i)M{A)) ^ 2(0(1)0-1(2)0(3)0-1(4)) + (2 ^ 3) + (3 ^ 4) (58) 

The permutation 2 <-> 3 maps z to 1/z, whereas 3 4 maps it to z/ {z— 1). The analytic continuation of the conformal 
blocks (|5^ ) to these regions can be performed directly by using hypergeometric function identities, or more easily by 
modifying the original screening charge integration contours appropriately. We find 

F^{l/z) = -z3/2Ff(z) 
(1/z) = -z^/^Fliz) 
/'2ni/2)--4z3/2F(l/2,5/2;2;z) 

F2^(l/z) = 2V2z=^/2^(l/2, 5/2; 2; 1 - z) (59) 

and 
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(^) = 2a/2(1 - z)3/2f(1/2, 5/2; 2; z) 
^i'(l^) = -4(1 - zf'^F{l/2, 5/2; 2; 1 - z) 



(60) 



These, together with ( |56| ) and (p9D, aUow us to write the exact expression for the full four-point M correlator in our 
model: 



(M(1)X(2)7W(3)M(4)) 



1 



where 



1-z 



-.■Hi 



1 - z 



^21 



1-Z 



22 



Huiz, z)= i/;2(z, z) - z|zr^i/;2(i 4) - (1 - z)|l - zp^i^^^l 



z 



,1 1, 



z- 1 z- 1' 

z z , 



H2l(z,z)=i/^l(z,z)-z|zr3^^l(-,^)-(l-z)|l-Zr3ff^l( 

zz z— Iz— 1 



7^22(2, 2-)= H',,iz,z) + \z\-^H',,{-, 4) + |1 - ^r^i/^2(^, ^ 

Z Z Z Z 1. 



(61) 



(62) 



with the correlators i/^^ given by 



HUz,z)= Ftiz)Fhz) + F^iz)Fnz) 



(63) 



From Eq. (^ij) we extract the OPE of M with itself, which in turn determines the scaling of products of local densities 
of states. 

VI. FUSION RULES AND OPERATOR ALGEBRA 

We will start here by studying the fusion rules of coset operators h{z, z). There are two cases of interest: the fusion 
hh^^, and hh. These two cases can be obtained from the coset correlator 



(/laibi {zi,Zi)h^,^l^ (22, Z2)/la3b3 (^3, 23)^^4(^4, Z4)) 



Zl2Z3iZl4Z23 



Z13Z24 



3/2 



5^/4- \F,^{z)F^{z) + a^b 



(64) 



^,3 



by respectively considering the limits zi — > Z2 in the above correlator, and in the one with 2 and 3 permuted. These 
correspond to taking z ^ in the appropriate solutions, which involves in the second case the analytic continuations 
z l/z. 

Using the formulas contained in the appendix, we find that the dominant contributions come from iJn, associated 
with the singlets (in the second case, this means contracting 1 with 3, 2 with 4). This yields (we here omit for 

simplicity the index structure) 



h{zi,Zi)h ^(Z2,Z2) ~ |2l2|^^^ 



^A{Z2) + ^A{Z2) +41+ 2^(Z2, Z2) + In |Z12|C(Z2, Z2) 
■212 212 



(65) 



in which X is the identity operator. The correlators of the operators appearing in ( pq ) are 

(A(zi)A(z2)) ^ zf,; {AizMiz2)) - ^2 

(i)(zi,Zi)l)(z2,Z2)) -ci - ln|zi2| 
(^(zi,Zi)C'(z2,Z2)) -1 
(C(Z1,Z1)C(Z2,Z2)) =0 



(66) 



(67) 
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where ci is some constant, unimportant for our purposes. 

Some comments are in order here. We notice that the most relevant operators appearing in the OPE (|6^), A{z) 
and A{z), possess conformal weights (—1,0) and (0,-1) respectively. Usually, the fusion rules for WZNW models 
would imply that the adjoint operator, whose conformal weights are 



N 



N- 



^ -1 



(68) 



should appear in the OPE (|65|). But the term in the four-point function pointing to the presence of such an operator, 
corresponding to a term inside the brackets of (65), does not appear since we are not allowed, by the 

requirement of single valuedness, to multiply the b contour solutions in the holomorphic and antiholomorphic sectors 
together. This requirement cannot be seen from the chiral conformal algebra: it comes out of the solution for the full 
correlator. 

Thus, single valuedness in the logarithmic conformal field theory that we are considering here changes the scaling 
behaviour expected from considerations not taking it into account, as in In fact, the operators ( ^6| ) exactly 
correspond to the local fields considered there. 

We can now very easily derive the OPE of OO^^. It has the exact same logarithmic form as above, with the 
exception that the proper fermionic contractions modify the weight in front: 



0(Z1,Z1)0-1(Z2,Z2) 



1 



This gives the new scaling formula 



4/3 

Zl2 



4/3 

Zl2 



A{z2) + 4X + 2D{z2, Z2) + In |zi2|C(z2, Z2] 



(69) 



([00-i](l)[00-i](2)) cxz2 + ^2 

to be contrasted with what was obtained in The other scalings mentioned there will also be modified for = 2. 
The correct expressions are given below. 

The OPE of the local DOS are extracted from Eg. (pi]) (we normalize the M two-point function). First, we can 
write the behaviour of (|6l|) as zi — Z2 = e ^ 0: 



[XX](2)7W(3)7W(4)) ^ 



1 



6^34 1 



1 + 7 In I e I + 7 In 



234 


+ ... 


2:23^24 





(71) 



in which 7 is a nonzero constant whose exact value is irrelevant for our purposes (it could be obtained by carefully 
expanding (|6l|)). The crucial fact is that the operator in the symmetric representation, with weights (—1,-1), does 
not appear here again, like the adjoint operator in (^5|). There is one quick way of seeing this last result. As z — > 0, 
~ —z. Thus, summing the conformal blocks for z and amounts to adding the OPE ( |69| ) for z and —z, 
which makes the antisymmetric 1/ z part disappear. The consequence of this, when one considers the point-splitting 
procedure leading to Al', is that negative- dimensional operators simply do not come about in the multifractality for 
SU{2). 

Postulating for M. an OPE of the form (we remove the tildes from the logarithmic operators, since we use a new 
normalization) 



M{l)M{2) 



then yields 



1 



\Zl2 



1/2 



X + i?(2) + iln|zi2|C(2) 



(C(1)M(2)7W(3)) = 27- 
{D{l)M{2)M{i))^^- 



1 



^23^/2 
1 



F23 



1/2 



In 



(i?(l)Z?(2)) =-27ln|zi2| 
{D{\)C{2)) - 27 
(C7(1)C(2)) = 



(72) 



2:23 



Z12Z1Z 



(73) 



11 



in turn leading to the OPEs 

i?(l)i^(2) = -27ln|zi2| + ..., 2?(1)C(2) = 27 + ..., C(1)C(2) = 

D{l)M{2) = --f\n\zi2\M{2) + C{1)M{2) = 2-fM{2) + ... (74) 

which wiU be used later to infer the behaviour of the multifractality. 



VII. iV > 2 CASE 

Let us now outline what happens for the N > 2 case. We have seen in the previous section that, for SU{2), the 
spectrum of operator dimensions is cut right at the beginning by the logarithmic locality conditions. We can wonder 
if the same holds in the N > 2 case. 

The conformal blocks, solutions to the SU{N)-2N Knizhnik-Zamolodchikov equations, are {F[{z) — [z{l — 

F^{z)^ V2N{1 - z)F{2-l/N,2 + l/N]2-z) 
Fl{z)^ (1 - z)F{2 - 1/7V, 2 + 1/N- 3; 1 - z) 
F^{z)= zF{2 -l/N,2 + 1/iV; 3; z) 

F^{z)= V2NzF{2 -l/N,2 + l/N;2;l- z) (75) 

For the case = 2, we had that the analytic continuations to 1/z were of logarithmic type. Here, however, the 
correct single- valued construction is like the standard one |3^ , since the analytic continuations 

F^{l/z)= -V2Nz-^{l - z)F{2 -l/N,2 + l/N; 2; 1/z) 
Fl{l/z)= -zi+i/^(l - z)F[l + 1/iV, 2 + 1/iV; 3; z) 
F^{l/z)^ z-^F{2 - l/N, 2 + 1/iV; 3; 1/z) 

f'|(l/z)= V2Nz'^+'^/^F{l/N, 2 + l/N;2;l- z) (76) 
are free of logarithms near z = for 7^ 2. On the other hand, we have for the blocks 



Ff (^)= V2iV(l - z)i+i/^F(l/iV, 2 + 1/iV; 2; z) 
P^ij^)= (1 - ^)"'^(2 - 2 + 1/N; 3; ^) 
^2 (i3t)= -(1 - zy+^/^'zFil + 1/N, 2 + 1/N; 3; z) 

^2(7^)^ ~V2N{1 - z)-hF{2 - 1/N, 2 + 1/N; 2; -^). 
^ 1 — z 



The full correlator reads in this case 
{M{1)M{2)M{3)M{4)) ■ 

where 



|l-z|^ 



'Hll + ="^12 + ^21 + T; F7^22 

1 — z 1 — z |1 — z|^ 



Hii{z,z)= H[^{z,z) + 


\z\- 


-2+^/^^i/li(i,i) 

z z 


+ ll 


Hl2(z,z)== H[2{z,z) - 


z\z 






n2i{z,z)= H'2^{z,z) - 


-z\z 


-4(l-l/iV^)^,^(l 

z 


z 


'H22{z,z)= H22{z,z) + 


\z\- 


-^^^-^/^')H'22{-. 

z 


z 


with the correlators H[j given by 









z 



'z- 1' z- 1-^ 

1 IAf^\ , , Z Z 



'"■^z-V z-V 
2^ z- 1' z- 1' 

r z-i' 



(77) 



(78) 



(79) 
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H^iz, -z) = F^{z)F^{z) + F^{z)F^{z). (80) 

As for the S'C/(2)_4 case, the b contour solutions have logarithmic behaviour near z = for z and -j^j (see appendix). 
Moreover, this logarithmic behaviour is exactly of the same form, involving the same powers of z. Thus, again, 
the operator in the adjoint representation, having conformal weights (—1,-1), does not appear in the OPE of the 
decoupling field with its inverse, because of single valuedness. 

However, when one considers the OPE of h with h, for which the analytic continuation to 1/z is necessary, the 
difference with the S'[/(2)_4 becomes clear. In fact, the analytic continuations of F to 1/z, in contrast to the N = 2 
case, are completely regular for > 2, as we have seen above. Expanding the full correlator ( |78| ) around z = 0, we 
find that the operators in the symmetric (S) and antisymmetric (A) representations do appear, as proposed in 
These have conformal weights 

Thus, we now have a negative dimensional operator on the physical level, in contrast to the N = 2 case. This operator 
has conformal weights 

h^h=^{l~2/N) (82) 

For N = 2, this vanishes, in accordance with the work of the previous sections. 

Logarithms in conformal field theories originate from degeneracies in the spectrum of conformal dimensions. When 
two or more operators have weights which become degenerate, they become distinct again by incorporating powers 
of Inz, which is a weight zero object. For the N = 2 case, three operators in the theory were of weight zero: the unit 
operator, the first descendant of the adjoint operator, and the operator in the antisymmetric representation. This high 
level of degeneracy has produced the logarithmic degeneracies which, through the requirements of single-valuedness 
(locality) of the correlator, have led us to discover the termination of the spectrum. 

In the case of A^ > 2, we have just seen that the termination does not appear immediately. We can however 
speculate about the moment at which the termination should occur. This relies on basic arguments of representation 
theory. For SU{N)^2N , the conformal weights of the operators in a given representation characterised by a Young 
tableau with boxes in the i-th row is given by the associated Casimir 



i=i 1=1 



as (considering holomorphic and anti-holomorphic) 

hh} = ^{/.} = (84) 

Now it is easy to check that for a completely antisymmetric representation of SU{N), i.e. fi = 1,« = 1,...,A^, 
the value of the Casimir vanishes. The associated operator thus has vanishing conformal weights, which makes it 
degenerate with the unit operator and the descendant of the adjoint operator, like in the N = 2 case extensively 
presented above. This operator in the completely antisymmetric representation is generated only by a product of A 
fundamental representations (our basic operator), so will appear only for q = N. We thus speculate that, for N > 2, 
the additional conformal weight degeneracy at q — N leads to logarithmic degeneracies which, exactly in the same 
way as previously, lead to the termination of further scaling. 



VIII. TERMINATION OF THE PARABOLIC MULTIFRACTAL SPECTRUM 

The relationship between the operator dimensions in the CFT and the multifractal spectrum were derived with the 
help of real-space renormalization group arguments in where all the definitions and notations used here are put 
forward. Three types of multifractal spectra are defined: T(q), which is the disorder averaged logarithm of the inverse 
participation ratio of normalized wavefunctions; r* (q) (the one we calculate here), which is the logarithm of the ratio 
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of disorder averaged unnormalized wavefunctions and disorder averaged normalization to the appropriate power, and 
T**(g), which is the logarithm of the disorder averaged participation ratio of normalized wavefunctions. There exist 
strong constraints on the first of these, T{q). It can only be a monotonously increasing function of q, which forces the 
cut of the parabolic approximation at a given point. The other functions have not been shown to have to obey to 
such conditions. 

The only multifractal spectrum that can be calculated with the help of the present theory is T*{q). It is related 
to the conformal weights (hq, hq) of the most relevant operator contained in the point-splitting definition of M"^ 
according to 



T*(q) = (q-1)(2-Ai) + A,-Ai 



(85) 



where = hq + hq. These weights are obtained by considering the decay of the correlator 



{1)^(2)) 



^12 ^1 



-12 



(86) 



Let us first recall the arguments used by ||] to obtain these weights. The operator O (^) is made from the product 
of free fermions and of a field in the fundamental representation of SU{N)-2N ■ The fusion of the fermions is 
straightforward, but the fusion of the WZNW primary fields is a bit less obvious. Building the Young tableaux 
associated to the tensor product of q fundamental representations and evaluating the quadratic Casimir eigenvalues 
c,n of all the resulting irreducible representations, one obtains a list of operators having conformal weights [^sP] 
hm = —Cm/N. It is then easy to identify the most relevant operator as the one in the completely symmetric 
representation, which has hs = —^^q{q + N). The resulting combined fermionic and WZNW weights lead to the 
provisional result M 



hq = hq 



N -I 
~2N^ 



[q^ + Nq) (provisional) 



and led to the corresponding scaling exponents t* (q) given by 



(87) 



r*(q) = (q-l)[2-^q 



(provisional) 



(88) 



One might question the fact that this t* (q) becomes arbitrarily negative for large values of q. It attains its maximum 
value at qmax = N(\ — + 1/2 > TV, after which the parabola bends back down. Even though it is in principle 

not impossible that T*(q) could take on negative values (see the discussion in [^), a limit to this behaviour would be 
most welcome. 

However, we can see from the previous sections that the straightforward WZNW fusion rules are not applicable 
in our case. Namely, we can perform for example the explicit point-splitting procedure leading to for the case 
N — 2 from the fusion rules (ItI): 



M^(z) = lim M(z + a)M(z)M(z - a) = lim \a\-^^^ 



I + D(z) + ^ln\a\C(z) + . 



(lim |a| 



-1/2 



)M(x) 



M(z-a) 
.. - M(x) 



(89) 



showing exphcitly that {M^(l)M^(2)) 



\Zl2 



i.e. that h-i ~ h-i — 1 /8, in contrast to the results of 



-3/8. Generalizing to arbitrary q (still for the specific case of = 2), we obtain that 



for which 



cx X, 



M 



2p+l 



CX M 



(90) 



as far as the correlators (A^'^(1)A^^(2)) are concerned. This implies the explicit termination of the parabolic law 
for the multifractality in the case of SU(2) vector potential randomness, which can be understood as a vanishing 
of OPE coefficients associated to operators with more relevant dimensions. This type of termination seems special 
to the SU(N) non-Abelian disorder case, since the Abelian multifractality, when computed from CFT, only involves 
expectation values of vertex operators of Gaussian distribution which do not carry the same conformal block structure 
as the ones that we have here. Note also that in the case of the gl(l, 1) WZNW model [^, fields with arbitrarily 
negative conformal dimensions can be obtained from the fusion of a sufficient number of fundamental representations. 
These interesting differences still need to be more properly understood. 
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Finally, from the results and discussion of the previous section, we can conjecture that the multifractal spectrum 
( ^ ) for general N should be complemented by a condition of the sort 

q < (91) 

where qc ^ N is the order at which the proper treatment of the locality conditions in the presence of logarithmic 
degeneracies kill off further scaling. Thus, above this critical value of q, the parabolic law will not hold anymore. 
From our explicit results for SU{2), and comparing with the results from M, we can say that 2 < qc{SU{2)) < 3. 



This is consistent with the exact solution for T(q) itself |40 



IX. WHAT IS WRONG WITH REPLICAS ? 

In general, the theoretical treatment of disordered systems is a notoriously difficult task in view of the fact that it 
is necessary to average the logarithm of the partition function over the statistical ensemble in order to obtain physical 
results. Two general approaches have been developed to perform such types of calculations. First and foremost, 
the supersymmetry (SUSY) approach is a mathematically sound technique that makes use of commuting and 
anticommuting fields simultaneously, making the partition function Z equal to one by definition. Its only setback is 
its limited applicability: interacting disordered systems cannot be represented by a SUSY field theory. 

The alternative technique is the well-known replica approach, introduced in |38| , which allows to rewrite the average 
of In Z as a tractable object by making use of the identity 

hTZ = lim ^LzJ. (92) 

r->0 r 

i.e. by considering a replicated theory with r species of the original field. One then hopes that the physical properties 
of the disordered system are faithfully represented by the analytical continuation to r ^ 0. However, the main problem 
is that the replica approach is not mathematically well-defined. Indeed, even though it has been used successfully 
quite often, the pathologies persist (for example, the systems obtained by choosing a positive or negative number of 
replicas are markedly different, making the limit r ^ completely ill-defined). It has been argued in js^ that failure 
of the replica trick may occur whenever the theory for a general integer value of r possesses different symmetries than 
the one at r = 0. 

We are thus now in a position to compare our present exact results with the previous ones obtained with the replica 
approach The conformal blocks obtained there were solutions to SU{0)n equations, and different to the ones 
obtained here. The OPE of the physical field (Q in [^) was given by 

M{l)M{2) ^ j;^^ [I - 4m + C(2) In \z\) - z(5(2) + C(2) In |z|) + ...] (93) 

where the chiral-like logarithmic operators had correlators given by the expressions 



{Dil)D{2))^^[ln\z,2\ + c] 

^12 

(Dil)C{2)) ^ i- 

(C(1)C(2)) = (94) 

with similar expressions for C, D given by z substituted by z. Although the basic two-point function generated from 
this OPE will be the same as the one obtained through above, the higher-point functions will show marked differences 
with the exact ones, which makes all further arguments unreliable. 

Thus, although the replica approach successfully accounts for the basic dimensionality of the operators, it does not 
correctly describe their higher-level correlations. The failure of the replica solution may be understood in terms of 
the interpretation put forward in [ p9| , since the essential ingredient in the replica solution was the SU{r) symmetry 
(associated to the r different flavours of replicated fermions), which obviously is a different symmetry at r = than 
at integer r. 



15 



X. CONCLUSION 



Let us give a brief summary of our results. 

(a) We have demonstrated explicitly that for the SU{2) group, the basic WZNW fusion rules are modified by the 
proper treatment of the locality conditions in the presence of so-called logarithmic operators, such that operators 
with negative conformal dimensions do not appear in OPE of the local density of states that we are interested in for 
multifractality purposes; 

(b) We have provided arguments based on the nature of the conformal field theory and its operator content, to 
conjecture that the mechanism for the termination of the multifractal spectrum that we propose should have occurred 
as well before q = N a,t the latest for general SU {N) . In these theories, there will thus be a finite string of physical 
operators with negative dimensions influencing the multifractality, but not an infinite set as previously believed. The 
termination for > 2 is conjectured to occur by the above mechanism before the scaling exponents r* {q) reach their 
maximum value, which provides evidence that the multifractal spectrum r* {q) obeys the same type of inequalities as 
the spectrum T(q) defined with normalized wavefunctions; 

(c) We have shown that the replica approach fails to reproduce the correct nonperturbative multipoint correlation 
functions in the conformal limit. 
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XII. APPENDIX 1: THE HAAR MEASURE 

For the sake of clarity, we here include the derivation of the Haar measure on SU'~'' (2)/ SU{2) as can be found in 

il- 

For a function f{g^) G Cc{G'~''), where Cc{G'~^) is the space of continuous functions of compact support on G*-^, 
there exists a function D{k, a, n) such that the group integration, after the Iwasawa decomposition, can be written as 

f{g'^)dg^ = f{kan)D{k,a,n)dkdadn (95) 

Jkan 

Since is unimodular (i.e. its left invariant measure is also right invariant), the LHS does not change when we 
replace fig^) by f{kiani),ki £ K,ni £ N. Then, it follows that D{k^^k, a,nn^^) = D{k,a,n), so D{k,a,n) is a 
function 5{a) of a only. Therefore, for ai e A, we have 

/ f{g'^)dg^ = / f{g^ai)dg = / f{kanai)S{a)dkdadn 
Jc^ JG<^ Jkan 

f[kaai{a^^nai))6{a)dkdadn — / f{ka{a^^nai))5{aa^^)dkdadn 
KAN Jkan 

f{kan)S{aai^)dkdaI{ai)*{dn) (96) 

KAN 

where /(ai) is the automorphism n oiria^^ of N, and I*uj denotes the transform (puUback) of uj by /. We have 
made use in the last step of the equality 

/ (/o$-i)cc>- / /$*^ (97) 
Jn Jm 

in which AI and N are two oriented manifolds, and $ is a diffeomorphism between them. 
Finally, by looking at 
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1 [/il+^/i2]p> 



(98) 



we see that I{ai)*{dn) — e^'^^ d^id^2 ~ e^"^M/i+(i/i_, and choosing ai = a produces the Haar measure (0). 



XIII. APPENDIX 2: ANALYTIC CONTINUATIONS 



We here give for reference the analytic continuation formulas used in the text, as well as the expansions of the 
relevant functions around z — Q. 

The second solutions for the conformal blocks, labeled by their contour 6, are hypergeometric functions of the 
argument 1 — z. When z ^ 0, these behave logarithmically. The proper expressions become 



TT 



4/3 ^ -FTW ^ _ 4/3 + (1 _ z)Ki^{z) 
z 2V2 



16/3 



+ F;^{z)\nz + Al-i + zKi2{z) 



F(l/2,5/2;2;l-z) = — 

TT 



-4/3 



F(l/2,5/2;2;z)lnz + ifoi(2) 



where 



(z)= £ [^(« + * + 1/2) + ^{n + 5/2) - ^.(n + j + 1) - + 1)] 



ri=0 



(n + j)!n! 



Near 2 = 0, we can then perform the expansions 



F^{z) = 2^/2[l + 7/8z + ©(z^)] 



4/3 



+ lnz-4/3 + is:ii(0) 



+ 0(zlnz) 



F2°(2) =2 + 5/4z2 + 0(23) 



V2 

TT 



16/3 



4/3 



O(zlnz) 



(99) 



(100) 



F(l/2, 5/2; 2; z) = 1 + 5/8z + 0{z^) 



F(l/2,5/2;2;l-z) = — 

TT 

used in the determination of the OPEs in the text. 
For the N > 2 case, we have for example 



-4/3 



+ lnz + /foi(0) 



0(z In z) 



r(i - i/iv)r(i + i/N) 

with J(z) some function regular as z ^ 0. 



Ar2 



(A^2 _ 



+ F{'(z)lnz- 



iV2 - 1 



+ (l-z)J(z) 



(101) 



(102) 
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